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1. Introduction 

Experiments at RHIC start with a baryon rich environment; hence they naturally have a non- 
zero chemical potential. The finite temperature field theory formalism easily admits a chemical 
potential; however, it leaves us with a complex action, and we can no longer use importance sam- 
pling. This results in the well known sign problem. If the chemical potential is small, we can 
employ the Taylor expansion method [[I]]. This method requires simulations only at zero chemical 
potential. Here we employ this method to study QCD with three dynamical quarks using the asqtad 
action [§] that we have already extensively studied at non-zero temperature, but without a chemical 
potential. Previously, we have performed a study with nonzero chemical potential only for N t =4 



2. Methodology 



We briefly review the formalism and methods that are detailed in Refs. [|T|] and ||[|. Physical 
quantities of interest are Taylor expanded in the chemical potentials (in physical units) fij and p.}, 
for light and strange quarks, respectively. We drop the bar when referring to the chemical potential 
in lattice units. For example: 



VT 3 



n,m=0 

Only terms with n + m even appear due to CP symmetry. 
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For the interaction measure, 
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Temperature dependendent coefficients c nm (T) and b nm (T) are combinations of observables that 
can be calculated on non-zero T ensembles, but with zero chemical potential. We Taylor expand up 
to sixth order. To compute all the required terms, 40 fermionic observables have to be determined 
using stochastic estimators, as well as several gluonic observables Ensembles are generated on 
a line of constant physics with mi = 0.lm s and m s approximately the physical strange quark mass. 
Our previous work used lattices with N t = 4. We now use N t = 6 and compare with the coarser 
lattices. Before we present our results, it is interesting to compare the free theory for different ./V, 
to see how the continuum limit is approached. (See Fig. [I] [gp.) 

Turning to the interacting theory we show the unmixed coefficients for the pressure in Fig. ||. 
There is considerable structure at low T and then an approach to the Stefan-Boltzmann (SB) limit 
above the cross-over temperature. Also, the higher order coefficients are small, but their errors 
grow rapidly. Note that the errors are better controlled for N t = 6 (red) than they were for N t = 4 
(black). 

In Fig. |3| we show the mixed coefficients for the pressure. Similar figures are available for the 
coefficients that are relevant for the interaction measure. Due to lack of space, we will not show 
them here. 



2 



QCD equation of state at non-zero chemical potential 



Steven Gottlieb 




.-. standard cjt^ SB 
■ ■ standard c 20 /c, SB 
^» standard c 4| /c WSB 

^ Naikc O(/ c oo,SB 
b — a Naik c 2Q /c 20 SB 
Naikc 4 „/c 40SB 
cm p4 %/%), SB 
n-s p4 c 20 /c 2ft SB 

c 4(/ c 40. SB 



N. 




Figure 1: Several expansion coefficients as a function of N t for the free theory [||]. On the left, for non- 
vanishing coefficients, we show the the ratio to the Stefan-Boltzmann limit. On the right, for C6o, we show 
the value itself. 
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Figure 2: Unmixed Taylor coefficients c„o and cq„ as a function of temperture. New results for N t = 6 are 
shown in red; black is used for N t = 4. 
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Figure 3: Mixed Taylor coefficients c mn as a function of temperture. New results for N t — 6 are shown in 
red, black is used for N t = 4. 



3. Results 

With the coefficients in hand, we can calculate interesting quantities, such as pressure, inter- 
action measure, energy density, quark number density, and quark number susceptibility. Due to 
non-zero c n \{T) terms a non-zero strange quark density is induced even with /I/, = 0. To study the 
n s = plasma, we must, therefore, tune /!/, as a function of /!/ and T. 

In Figs. |] and |5| we show how the pressure and interaction measure change as a function of T 
for selected values of \ljJ '. In each figure, we display the zero chemical potential case on the left. 
We compare results for N t = 6 with our prior results for N t = 4. We find that the change in pressure 
is somewhat smaller compared with our previous results. For the interaction measure, the errors 
are fairly large, but there also seems to be a reduction there. Figure |6| shows the energy density and 
change in energy density due to chemical potential. In this case, the differences between N t =4 and 
6 are small. In Fig. [7|, we show the light quark number density and the quark number susceptibility. 
We note that both the number density and susceptibility are somewhat smaller with N t = 6 than 
they were for N t = 4. 

One particularly interesting quantity is the isentropic equation of state (EOS). In a heavy-ion 
collision, after thermalization, the system expands and cools with constant entropy. Thus, we would 
like to find the EOS with fixed ratio of entropy to baryon number. The appropriate ratio of s/ns 
for AGS, SPS and RHIC are 30, 45 and 300, respectively. To carry out this calculation we must 
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Figure 4: The pressure with zero chemical potential (left) and the change in pressure due to chemical 
potential (right). 




Figure 5: The interaction measure with zero chemical potential (left) and the change in interaction measure 
due to chemical potential (right). 




Figure 6: The energy density with zero chemical potential (left) and the change in energy density due to 
chemical potential (right). 
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Figure 8: Isentropic pressure (left) and interaction measure (right) for selected values of S/ns appropriate 
to AGS, SPS and RHIC. 

find trajectories in the {pii,}Xh,T) space with n s = and s/n B as stated above. In Fig. |8| we show the 
isentropic pressure and interaction measure. We also show the isentropic energy density, and quark 
number density in Fig. [| and the isentropic light and strange quark susceptibilities in Fig. [l0| 

4. Conclusions 

We have extended our sixth order Taylor expansion study of thermodynamics with chemical 
potential toward the continuum limit by going from N t = 4 to 6. After computing the expansion 
coefficients relevant for both pressure p and interaction measure / we can compute a number of 
interesting quantities. We observe modest lattice spacing effects, with the quark densities and 
susceptibilities, and the effect of chemical potential, smaller at the smaller lattice spacing. In 
addition, we have calculated the isentropic equation of state, which is particularly relevant for the 
phenomenology of relativistic heavy-ion colliders. For both values of N t we find rather smooth 
behavior for the isentropic variables indicating that experiments are far from any critical point in 
the pi-T plane. 

It would be interesting to extend this work to yet smaller lattice spacing and to go to lighter 
quark mass. 
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Figure 9: Isentropic energy density (left) and light quark number density (right) for selected values of S/ns 
appropriate to AGS, SPS and RHIC. 
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Figure 10: Isentropic light quark number susceptibility (left) and strange quark number susceptibility (right) 
for selected values of S/ns appropriate to AGS, SPS and RHIC. 
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